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Abstract. Numerous investigations on the fundamental properties of low-mass stars using 
eclipsing binaries indicate a strong discrepancy between theory and observations that is still 
awaiting explanation. Current models seem to predict radii for stars less massive than the 
Sun that are some 10% smaller than observed, while their effective temperatures are some 
5% larger. Here we discuss recent new observational data that are relevant to this issue and 
review the progress made in understanding the origin of the important differences with the- 
oretical calculations. Notably, we provide evidence that stellar activity may be responsible 
for the mismatch between observations and theory through two different channels: inhibi- 
tion of convection or effects of a significant starspot coverage. The activity hypothesis is put 
to a test with observational diagnostics and some of the consequences of the large starspot 
coverage are evaluated. We conclude that stellar activity likely plays a key role in defining 
the properties of active low-mass stars and that this should be properly taken into account 
when investigating young, active stars in clusters or star-forming regions. 
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1. Introduction 

Evidence collected over the past years con- 
vincingly shows that the components of low- 
mass eclipsing binary stars are not adequately 
reproduced by current evolution models. The 
observations yield stars that are systematically 
larger and cooler than theoretical calculations 
by about 10% and 5%, respectively. But, in 
contrast, these binary components are found 
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to have luminosities that are in good agree- 
ment with those of single stars and also with 
model predictions. This is the reason why stud- 
ies focused solely on the mass-luminosity re- 
lationship did not reveal any mismatch be- 
tween observations and theory and, further, 
obtained very tight mass and luminosity cor- 
respondence when using infrared photometry 
(e.g.. lDelfosse et al., 2000) . 

A revision of the available observational 
evidence and some hints about the possible res- 
olutions of the discrepancies were provided by 
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Fig. 1. Left: M~R plot for low-mass eclipsing binary stars with empirical determinations. Right: Same but 
for stars with uncertainties below 3%. The solid line represents a theoretical isochrone of 1 Gyr calculated 
with the lSaraffe et al.. C1998,) models. 



iRibasI (|2006|) . We refer to that study for back- 
ground information and detailed references. 
Here we shall only discuss new results ap- 
peared since 2006 and also the progress made 
in understanding the origin of the difference 
between model predictions and observations of 
fundamental properties of low-mass stars. 

2. Recent observational data 

New low-mass eclipsing binaries have re- 
sulted in the past two years from dedicated 
monitoring of carefully selected targets 
from photometric databases of variabil- 
ity surveys. These have bee n pii blished 



d2006h and 
In addition. 



jLopez-Morales & Shawl 
iLo p ez-Moral es et alj (l2006 l). 
the low-rn as s ecli psing binary reported by 
iHebb et al.1 (l2006l) resulted from a deep 
targeted search in several open clusters. In 
other cases, low-mass eclipsing binaries were 
serendipitously discovered over the course 
of photometric monitoring campaigns with 
different scientific aims, and later analyzed 
specifically in detail, such as the obje cts stud- 
ied by Bavless & Orosz (2006), Yo ung etal] 
toOdi . and lBlake et alJ (l2007h . A mass-radius 
plot of all presently known low-mass stars in 
detached eclipsing binaries is provided in Fig. 
[T](left), while Fig. [1] (right) shows only those 
objects that have reported error bars in both 



masses and radii below 3%. The systematic 
offset of 5-10% between the observations 
and the 1 Gyr isoc hrone from the models of 
iBaraffe etal J (119981) is apparent. 



In addition to the "classical" eclipsing bi- 
naries, there has been an increasing number of 
discoveries resulting from follow-up of plane- 
tary transit candidates. In some instances, the 
object responsible for the transit was found not 
to be a planet but an M-dwarf secondary to a 
F-G-type primar y star. This i s the c ase of the 
recent study by Beattv etakl (|2007|) of one of 
the HAT network planetary candidates. These 
objects are single-line and single-eclipse bina- 
ries that directly provide a value of the stel- 
lar density but that require certain assumptions 
(e.g., orbital synchronization) to derive the ac- 
tual masse s and r a dii. A similar technique ex- 
ploited by iTorresI (l2007l) has provided masses 
and radii (dependent on the assumed effective 
temperature) for the M-type star GJ 436, which 
hosts a transiting exoplanet. A mass-radius plot 
of all the objects resulting from single-line ra- 
dial velocities and single-eclipse light curves is 
shown in Fig.|2l The only object in this sample 
with an error bar below 3% (making it a re- 
liable test of models) is GJ 436. Interestingly, 
it also shows the 10% radius differential with 
model predictions tToiTe&.2007,) . 
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Fig. 2. Same as Fig. [T] (left) but for single-line and 
single-eclipse systems. 



Finally, there is an additional source of fun- 
damental properties of low-mass stars from di- 
rect interferometric measurements. This was 
the method applied in the recent study of 
iBer^er et all (l2006l) . where the authors re- 
port measurements for six M dwarfs with the 
CHARA array and suggested a correlation be- 
tween radius differences and stellar metallicity. 

In this context it is also worth reviewing 
the progress made in understanding the prop- 
erties of the only known eclipsing binary with 
brown dwarf com ponents. Following the dis- 



covery analysis of Sta ssun et al. 



(120061) . the re- 
20071) has fur- 



cent analysis of Stass un et al.l 
ther refined the properties of the binary compo- 
nents and strengthened the case for a reversed 
temperature ratio, which is difficult to recon- 
cile with model calculations. 

3. Stellar activity hypothesis: a 
theoretical framework 



As discussed by R ibasI (l2006h and lTorres et alJ 
( i2006.) there is a property that distinguishes 
low-mass stars in close binaries from those 
that are single objects, and this is the presence 
of tidal interactions that force the component 
stars to rotate in orbital synchrony. Because 
of strong observational biases, all low-mass 
eclipsing binaries found so far have periods 
well below 10 days, and, thus, rotation peri- 
ods also below 10 days. Using the analysis of. 



e.g.. |Pizzolato et all (l2003l) it is easy to realize 
that M stars with such short rotation periods 
will experience high levels of magnetic activ- 
ity. This is also confirmed by the fact that most 
low-mass eclipsing binaries are strong X-ray 
sources and also they show intense emission in 
the Ha Balmer line. 

From these observational facts it is sen- 
sible to consider a hypothesis by which the 
differences between model predictions and 
eclipsing binary observations arise from the 
high levels of magnetic activity of the com- 
ponent stars. Careful consideration of the ef- 
fects of magnetic fields in low-mass stel- 
lar evolution models was carried out by 
Mullan & MacDonald (2001), using some pre- 
scriptions from Gou gh & Tavleil (119661) . The 
authors pointed out that the inclusion of mag- 
netic fields potentially has a moderate impact 
on the overall stellar properties. More recently, 
Chabrier et al. (2007) have performed a more 
realistic treatment targeted specifically to re- 
solving the current differences between obser- 
vations and models. The authors have consid- 
ered two scenarios. Firstly, a scenario that con- 
siders the effect of magnetic fields in caus- 
ing inhibition of the convective energy trans- 
port. And secondly, a scenario similar to that 
proposed by Lopez-Morales & Ribas (2005) 
and based on simple energy conservation ar- 
guments in a spot-covered surface. 

The study of IChabrier et aP (l2007h con- 
cludes that both mechanisms alter the proper- 
ties of the star sufficiently to explain the ob- 
served radius and temperature discrepancies. 
In the case of the inhibition of convection, this 
was tested by setting the mixing-length param- 
eter to lower values than that yielded by the 
standard solar model. The result on the struc- 
ture of the star is that it shows a higher radius 
and slightly lower effective temperature. Tests 
show that good agreement with the observa- 
tions is obtained for a mixing-length param- 
eter of 0.5 pressure scale heights in the case 
of the more massive stars of the sample (i.e., 
M > 0.6 Mq). However, less massive stars 
are little affected by the lower mixing-length 
parameter and agreement with observations is 
only obtained when decreasing to (possibly un- 
physical) values of 0.1 pressure scale heights 




Fig. 3. Comparison of the best known eclipsing binary systems with models using p values of and 0.3. 
Left: M~R plot. Right: M-T^g plot. 



or less. This i s well illustrated in figure 1 of 
IChabrieretaT]d20Q7h . 

In the case of the direct effect of starspot 
coverage, the analysis can be carried out by as- 
suming a new luminosity L' expressed as 



L' = (1 -p)AnR'^o-T'^ 



eff 



(1) 



where R' is the modified radius, T'^^ is the mod- 
ified effective temperature andyS is the factor by 
which starspots block the outgoing luminos- 
ity because of their lower temperature. From 
the calculations of mo dels with this modi- 
fied boundary condition. IChabrier et alJ (l2007h 
found that the radii of low-mass stars in eclips- 
ing binaries in the entire mass domain can 
be reproduced with a fi parameter of about 
0.3. Assuming a spot temperature contrast of 
about 15% (or 500 K), the results indicate that 
starspots cover approximately half of the stel- 
lar surface. In this scenario, the modified stellar 
effective temperatures are also found to agree 
with the observations as the total stellar lumi- 
nosity is nearly invariant. In Fig. [3] we show 
a mass-radius plot and a radius-effective tem- 
perature plot that illustrate the good agreement 
between observations and models when con- 
sidering the effects of starspots. 

Thus, although conclusive evidence on 
which of the two scenarios (or perhaps a com- 
bination) is the most reliable to explain the 
observations is still lacking, the results of the 
analysis strongly suggest that stellai" activity is 



a key element in understanding the properties 
of low-mass stars such as their radii and effec- 
tive temperatures. 

4. Stellar activity effects: test with 
observations 

In this section we consider various observa- 
tional evidence in the context of the stel- 
lar activity hypothesis. One of the straightfor- 
ward questions to address is the universality 
of the hypothesis. In other words, if activity is 
thought to be responsible for the observed ra- 
dius and effective temperature discrepancy in 
eclipsing binaries, do single active stars also 
show the same effect? Indeed, active stars have 
long been recognized to define a slightly off- 
set s equence in t he color-magnitude diagram 
(e.g.. IStauffer & Hartmann 1986). The results 
from the eclipsing binaries and the theoreti- 
cal studies discussed above allow for a fresh 
look at the problem. From the available evi- 
dence it can be assumed that the luminosities 
of low-mass stars are not signifi cantly affected 
by stellar acti vity. The studies bv[ 



(Ithis vol ume*) and Morale s et al 



Morales et alJ 
(l2007l) show 



that a comparison of active and inactive stars 
of the same luminosity indeed reveals a sys- 
tematic temperature offset. The effective tem- 
perature differential can be translated into a ra- 
dius difference that is of the same order as that 
found in eclipsing binaries. This result gener- 
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alizes the activity effects on stellar properties 
to any star, either single or binary. Full details 
on the analysis can be found in the references 
above. 

The ac tivity hypothesi s was also inves- 
tigated by iLopez-MoralesI ( l2007l) . who col- 
lected values of rotational velocities and, even- 
tually. X-ray luminosities of a sample of low- 
mass eclipsing binaries. Then, the radius dis- 
crepancies were searched for correlations with 
such X-ray luminosities, which are proxies to 
the overall stellar activity. The author identi- 
fied significant correlations in the binary sam- 
ple but not so in single stars with interfer- 
ometric radius measurements. However, we 
must caution the reader that the dynamic range 
of the X-ray lu minosities of single stars in 
iLopez-MoralesI ( 120071) is very small and the 
overall scatter of the X-ray values will mask 
any correlation. Thus, the lack of inconclusive 
evidence in the single star analysis is not sur- 
prising. 

There is yet another observational piece of 
evidence that adds to the discussion on the ef- 
fects of activity on stellar properties. This is the 
detailed analysis of the only known eclipsin g 
brown dwarf reported bv lStassun et alj (l2007h . 
The surprising conclusion of the analysis was 
a statistically significant temperature reversal 
(i.e., the more massive brown dwarf is also 
the cooler of the pair), which is not compati- 
ble with the pred ictions of models. H owever, 
a recent study by Reiners et al.l (12007 ') indeed 
shows that the more massive component is 
significantly more active than its less massive 
counterpart. Put in the context of our results 
above, it seems clear that higher activity would 
result in cooler temperature thus explaining the 
reversal found in the light curve analysis. 

But there are also results that seem to be at 
odds with the hypothesis of stellar activity as 
being responsible for the ~ 10% differential in 
stellar radii. The recent study by lTorresI (l2007h 
used the particularities of planetary transits to 
determine the physical properties of the host 
star GJ 436. This star is relatively inactive, 
with a value logLx/i'boi = -5.1 (compared 
with -3, which is typical of eclipsing bina- 
ries). Surprisingly, the author also found a 10% 
radius differential when compared with mod- 



els, in agreement with the results from eclips- 
ing binaries. More examples of such objects 
with transiting planets are desirable to draw 
any firm conclusions. However, in the mean- 
time, we call attention to the fact that the anal- 
ysis of GJ 436 is heavily based on an assumed 
stellar effective temperature, which is known 
to be poorly established in the low end of the 
main sequence, and thus potentially subject to 
large systematic errors. 

5. Starspot coverage 

One of the consequences of the starspot sce- 
nario discussed above is the need for a fi 
value of 0.3 and thus about 50% surface cov- 
erage of spots cooler than the photosphere by 
15% to explain the observed radius discrepan- 
cies. It is worth reviewing now if such high 
coverage is compatible with the observations. 
Interestingly, large spot coverages are intu- 
itively associated with large photometric varia- 
tions. This is indeed not the case. Photometric 
variations are only sensitive to the contrast be- 
tween different areas on the surface of the star. 
Thus, a heavily spotted star with a homoge- 
neous spot distribution would have its overall 
light level severely diminished but display no 
significant variations along the rotation phase. 

We carried out various simulations to in- 
vestigate if the eclipsing binary data in the 
form of light curves are compatible with the 
inferred surface spot coverage of about 50%. 
This was done by assuming different latitude 
distribution models and then evaluate the peak- 
to-peak magnitude variations. A detailed sum- 
mary of our investigations will be provided 
elsewhere but here we point out that the simu- 
lated light variations out of the eclipses and the 
depth of the ecUpses themselves are compat- 
ible with the observed light curves of several 
analyzed eclipsing binaries. An illustration of 
this can be seen in Fig. H for the case of the 
eclipsing binary YY Gem. 

6. Other scenarios 

Stellar activity is not the only scenario advo- 
cated to explain the reported differences be- 
tween observed and model calculated stellar 
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Fig. 4. Four top panels: Plots depicting two different latitude distributions of starspots covering 50% of the 
stellar surface. Note the different out-of-eclipse light drop and the resulting phase variations. Bottom panel: 
Observed light curve of the eclipsing system YY Gem in the R band. 



radii. iBerger et al.l(l2006l) analyzed a sample of 
stars with interferometric radii and found quite 
a strong correlation between radius differences 
and stellar metallicity. From this evidence the 
authors conclude that the mismatch between 
observation and theory could be due to a miss- 
ing source of opacity in the model calcula- 
tions. However, a reana lysis of an extended in- 
terferometric sample by lLopez-MoralesI ( l2007h 
failed to identify such strong correlation. In 
ad dition, the interf e romet ric radii determined 
by ISegransan et al.l (l2003h with VLTI seem to 
agree well with model predictions. The possi- 
ble relationship between radius differences and 



metallicity should be further investigated, both 
with additional data and with a statistically 
sound approach to evaluate the significance of 
the correlations found. 



7. Conclusions 

Without excluding at this point any other sce- 
nario, it seems clear from the tests discussed 
above that stellar activity plays a key role in 
defining the structure and radiative properties 
of low-mass stars. This could be very impor- 
tant in the context of young, and therefore 
active, low-mass stars, for example in clus- 
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ters or star-forming regions. Failure to account 
for the effects of stellar activity can lead to 
strong potential biases in the determination of 
the ages of these objec ts and th eir ensembles. 
As discussed by MoralesietalJ |2007), a sim- 
ple calculation using the temperature differen- 
tials found for single active stars indicates that 
ages of young clusters determined from (ac- 
tive) low-mass stars could be systematically 
underestimated by about 40%. Interestingly, 
this difference agrees well with the discrep- 
ancy between color-magnitude diagram ages 
and Li de pletion boundary ages in young c lus- 
ters (e.g.. lBarrado y Navascues et alll2004i) . 

Acknowledgements. The authors acknowledge 
support from the Spanish Ministerio de Educacion 
y Ciencia through the program for Acciones 
Integradas HF2005-0249 and PNAyA grants 
AYA2006-15623-C02-01 and AYA2006- 15623- 
C02-02, and from the French Picasso program 
11412SB. 



References 

Baraffe, I., Chabrier, G., Allard, R, Hauschildt, 

P. H. 1998, A&A, 337, 403 
Barrado y Navascues, D., Stauffer, J. R., & 

Jayawardhana, R. 2004, ApJ, 614, 386 
Bayless, A. J., & Orosz, J. A. 2006, ApJ, 651, 

1155 

Blake, C. H., & Torres, G., & Bloom, J. S., 
& Gaudi, B. S. 2007, ApJ submitted (astro- 
ph/0707.3604) 

Beatty, T. G., & Fernandez, J. M., & Latham, 
D. W., & Bakos, G. A., & Kovacs, G., & 
Noyes, R. W., & Stefanik, R. R, & Torres, 

G. , & Everett, M. E., & Hergenrother, C. W. 
2007, ApJ, 663, 573 

Berger, D. H., & Gies, D. R., & McAhster, 

H. A., & Brummelaar, T. A. t., & Henry, 
T. J., & Sturmann, J., & Sturmann, L., 
& Turner, N. H., & Ridgway, S. T, & 
Aufdenberg, J. R, & Merand, A. 2006, ApJ, 
644, 475 

Chabrier, G., Gallardo, J., & Baraffe, L 2007, 
A&A, 472, 17 

Delfosse, X., & Forveille, T., &Segransan, D., 
& Beuzit, J.-L., & Udry, S., & Perrier, C, & 
Mayor, M. 2000, A&A, 364, 217 



Gough, D. O., & Tayler, R. J. 1966, MNRAS, 
133,85 

Hebb, L., & Wyse, R. F G., & Gilmore, G.,& 
Holtzman, J. 2006, AJ, 131, 555 

Lopez-Morales, M. 2007, ApJ, 660, 732 

Lopez-Morales, M., Orosz, J. A., Shaw, J. S., 
Havelka, L., Arevalo, M. J., Mclntyre, T, 
& Lazaro, C. 2006, ApJ, submitted (astro- 
ph/0610225) 

Lopez-Morales, M., Shaw, J. S. 2006, in 
7th Pacific Rim Conference on Stellar 
Astrophysics, ASP Conference Series, eds. 
Kang, Y. W, & Lee, H. W., & Cheng, K. S., 
& Leung, K. C. (astro-ph/0603748) 

Lopez-Morales, M., & Ribas, I. 2005, ApJ, 
631, 1120 

Morales, J. C, & Ribas, I., & Jordi, C. in XXI 
Century Challenges for Stellar Evolution, 
Memorie della Societa Astronomica 
Italiana, Eds. S. Cassisi, & M. Salaris (this 
volume) 

Morales, J. C, & Ribas, I., & Jordi, C. 2007, 
A&A, in press (arXiv:071 1.3523) 

Mullan, D. J., & MacDonald, J. 2001, ApJ, 
559, 353 

Pizzolato, N., Maggio, A., Micela, G., 

Sciortino, S., & Ventura, P 2003, A&A, 

397, 147 
Ribas, I. 2006, Ap&SS, 304, 89 
Reiners, A., & Seifahrt, A., & Stassun, K. G., 

& Melo, C, & Mathieu, R. D. 2007, ApJ, in 

press (astro-ph/071 1.0536) 
Segransan, D., Kervella, P., Forveille, T., & 

Queloz, D. 2003, A&A, 397, L5 
Stauffer, J. R., & Hartmann, L. W. 1986, ApJS, 

61,531 

Stassun, K. G., & Mathieu, R. D., & Valenti, 

J. A. 2006, Nature, 440,311 
Stassun, K. G., & Mathieu, R. D., & Valenti, 

J. A. 2007, ApJ, 664, 1154 
Torres, G. 2007, ApJ, in press (astro- 

ph/0710.4883) 
Torres, G., Lacy, C. H., Marschall, L. A., 

Sheets, H. A., & Mader, J. A. 2006, ApJ, 

640, 1018 

Young, T. B., Hidas, M. G., Webb, J. K., 
Ashley, M. C. B., Christiansen, J. L., 
Derekas, A., & Nutto, C. 2006, MNRAS, 
370, 1529 



